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Charm-Production in e+e Annihilation Around 4 GeV 

B. W. Lang 

University of Minnesota, 116 Church St. S.E., Minneapolis, MN 55455 

Using the CLEO-c detector at the Cornell Electron Storage Ring, we have measured inclusive and exclusive 
cross sections for the production of -D"*", and mesons in e+e" annihilations at thirteen center-of-mass 
energies between 3.97 and 4.26 GeV. Exclusive cross sections are presented for final states consisting of two 
charmed mesons (DD, D*D, D* D* , D'^DJ , D^'^ DJ , and Ds+D^") and for processes in which the charmed 
meson pair is accompanied by a pion. 



1. Introduction 

Hadron production in electron-positron annihila- 
tions just above cc threshold has been a subject of 
mystery and little intensive study for more than three 
decades since the discovery of charm. Recent devel- 
opments, like the observation of the y(4260) reported 
by the BaBar collaboration [i] and subsequently con- 

; firmed by CLEO-c [2] and Belle [l], underscore our 
incomplete understanding and demonstrate the poten- 
tial for discovery of new states, such as hybrids and 
glueballs. It is also clear that precise measurements 
of charm-meson properties are essential for higher- 
energy investigations of ^-flavored particles and new 
states that might decay into b. They also offer unique 

' opportunities to test the validity and guide the devel- 
opment of theoretical tools, like Lattice QCD, that are 
needed to interpret measurements of the CKM quark- 
mixing parameters Q. Any comprehensive program 
of precise charm-decay measurements demands a de- 
tailed understanding of charm production. 

Past studies of hadron production in the charm- 

\ threshold region have been dominated by measure- 
ments of the cross-section ratio R — a{e^e~ 
hadrons, s)/a{e^e~ ^.^ , s) that have been made 

over this energy range by many experiments Re- 
cent measurements with the Beijing Spectrometer 
(BES) 0] near charm threshold are especially note- 
worthy. There is a rich structure in this energy region, 

, reflecting the production of cc resonances and the 
crossing of thresholds for specific charm-meson final 
states. Interesting features in the hadronic cross sec- 
tion between 3.9 and 4.2 GeV include a large enhance- 
ment at the threshold for D*D* production (~ 4 GeV) 
and a fairly large plateau that begins at D^'^Dj 
threshold. While there is considerable theoretical in- 
terest 0, H, , there has been little experimental 
information about the composition of these enhance- 
ments. 

In this paper we describe measurements of charm- 
meson production in e"'"e~ annihilations at thirteen 
center-of-mass energies between 3970 and 4260 MeV. 
These studies were carried out with the CLEO-c de- 
tector at the Cornell Electron Storage Ring (CESR) 
[m in 2005-6. (Throughout this paper charge- 
conjugate modes are implied.) The principal objective 



of the CLEO-c energy scan was to determine the op- 
timal running point for studies of D+-meson decays. 
The same data sample has been used to confirm the 
direct production of F(4260) in e+e^ annihilations 
and to demonstrate y(4260) decays to final states in 
addition to tt+tt^J/V' [4I . Specific results presented 
in this paper include cross-section measurements for 
exclusive final states with Z?+ , and mesons and 
inclusive measurements of the total charm-production 
cross section and R. 



2. Data Sample and Detector 

The data sample for this analysis was collected with 
the CLEO-c detector. Both the fast-feedback analy- 
sis carried out as data were collected and the detailed 
analysis reported here are extensions of techniques de- 
veloped for charm-meson studies at the -0(3770) [l^ . 

An initial energy scan, conducted during August- 
October, 2005, consisted of twelve energy points be- 
tween 3970 and 4260 MeV, with a total integrated lu- 
minosity of 60.0 pb~^. The scan was designed to pro- 
vide cross-section measurements at each energy for all 
accessible final states consisting of a pair of charmed 
mesons. At the highest energy point these include 
DD, D*D, D*D\ D+D-, D*+ D' , and D*+ Df , 
where the first three contain both charged and neu- 
tral states. A follow-up run beginning early in 2006 
provided a larger sample of 178.9 pb~^ at 4170 MeV 
that proved essential in understanding the composi- 
tion of charm production throughout this energy re- 
gion. The center-of-mass energies and integrated lu- 
minosities for the thirteen subsamples are listed in 
Table m Integrated luminosity is determined by mea- 
suring the processes e+e~ e+e~, n~^^J-~, and 77, 
which are used since their cross sections are precisely 
determined by QED. Each of the three final states 
relies on different components of the detector, with 
different systematic effects. The three individual re- 
sults are combined using a weighted average to obtain 
the luminosity needed for this analysis. 

CLEO-c is a general-purpose magnetic spectrome- 
ter with most components inherited from the CLEO 
III detector [l^, which was designed primarily to 
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Table I Center-of-mass energies and integrated luminosity 
totals for all data samples for the CLEO-c energy scan. 



iicm (MeV) 


J L at (pb j 


3970 


3.85 


3990 


3.36 


4010 


5.63 


4015 


1.47 


4030 


3.01 


4060 


3.29 


4120 


2.76 


4140 


4.87 


4160 


10.16 


4170 


178.89 


4180 


5.67 


4200 


2.81 


4260 


13.11 



study B decays at the T{AS). Its cylindrical charged- 
particle tracking system covers 93% of the full 47r 
solid angle and consists of a six-layer all-stereo in- 
ner drift chamber and a 47-layer main drift cham- 
ber. These chambers are coaxial with a supercon- 
ducting solenoid that provides a uniform 1.0-Tesla 
magnetic field throughout the volume occupied by 
all active detector components used for this analysis. 
Charged particles are required to satisfy criteria en- 
suring successful fits and vertices consistent with the 
e'^e~ collision point. The resulting momentum res- 
olution is ^ 0.6% at 1 GeV/c. Oppositely-charged 
and vertex-constrained pairs of tracks are identified 
as K'g — > TT+TT^ candidates if their invariant mass is 
within 4.5 standard deviations (ct) of the known mass 
(~ 12 MeV/c2). 

The main drift chamber also provides dE/dx mea- 
surements for charged-hadron identification, comple- 
mented by a Ring-Imaging Cherenkov (RICH) detec- 
tor covering 80% of An. The overall efficiency for pion 
or kaon identification is greater than 90%, and the 
misidentification probability is less than 5%. 

An electromagnetic calorimeter consisting of 7784 
CsI(Tl) crystals provides electron identification and 
neutral detection over 95% of Att, with photon-energy 
resolution of 2.2% at 1 GeV and 5% at 100 MeV. 
We select tt" and t] candidates from pairs of photons 
with invariant masses within 3a of the known values 
[|(~ 6 MeV/c^ for tt" and - 12 MeV/c^ for rj). 

3. Event-Selection Procedures 

The procedures and specific criteria for the selec- 
tion of D^, and mesons closely follow previ- 
ous CLEO-c analyzes and are described in Refs. 
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Figure 1: A_B vs. Mbc in a Monte Carlo simulation of ex- 
pected final states at a center-of-mass energy of 4160 MeV, 
showing clear separation among the expected two-charm- 
meson final states. 



and [14|. Candidates are identified based on their in- 
variant masses and total energies, with selection cri- 
teria optimized on a mode-by-mode basis. We use 
only the cleanest final states for {K~n^) and 

(_ftr~7r+7r+) selection, since these provide sufficient 
statistics for precise cross-section determinations. For 

we optimize for efficiency by selecting eight de- 
cay modes: (/'tt"'", K*^K^, V'^'^ ^ VP^^ V"""*"; v'p'^^ 
ipp^, and KgK~^. Accepted intermediate-particle de- 
cay modes (mass cuts) are (f> K'^K~ (±10 MeV), 
X*o _^ x-TT- (±75 MeV), r]' t^tt+t:- (±10 MeV), 
and p+ -f 7r+7r° (±150 MeV). 

To determine the production yields and cross sec- 
tions for the final states accessible at a particular 
center-of-mass energy, we classify events based on the 
energy of a Di^g-^ candidate ( Ai? = £^beam — ) ^.nd 
its momen tum in the form of beam-constrained mass 
(^4c = \/Sbcam-|^i5(.,P)- Figure [I] shows the ex- 
pected behavior in a two-dimensional plot of Ai? vs. 
Afbc for a Monte Carlo simulation of CLEO-c data 
at 4160 MeV with about ten times the statistics of 
our data sample at that energy. There is clear sep- 
aration of events into the expected final states con- 
sisting of two charmed mesons. This separation was 
exploited during the scan run for a fast-feedback " cut- 
and-count" determination of event yields. It is also 
evident in plots of the momenta of charm-meson can- 
didates selected by cutting on candidate mass that 
the composition of final states can be analyzed by fit- 
ting the momentum spectra of D+ and D+ can- 
didates. Figure [2] illustrates this with the momentum 
spectra for K^n^ candidates within 15 MeV of 

the nominal mass both in the Monte Carlo sample of 
Fig.[T]and in 10.16 pb'^ of CLEO-c data at 4160 MeV. 
While no background corrections have been applied to 
these distributions, the structure of distinct Doppler- 
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Figure 2: Momentum spectra at 4160 MeV for D" 
K~'k'^ candidates with invariant masses within 15 MeV 
of the nominal value for Monte Carlo (top) and data (bot- 
tom). As described in the text, the concentrations of en- 
tries correspond to the dominant expected final states with 
two charmed mesons. 



smeared peaks corresponding to different final states is 
evident. The Monte Carlo and data show good quali- 
tative agreement, with concentrations of events corre- 
sponding to the dominant final states near 0.95 GeV/c 
(DD), 0.73 GeV/c {D*D) and 0.5 GeV/c {D*D*). 
The cross sections for all contributing final states can 
be determined by correcting the raw measured mo- 
mentum spectra like Fig.[2]for combinatoric and other 
backgrounds and then fitting to Monte Carlo predic- 
tions of the spectra. To achieve good fits, all signifi- 
cant production mechanisms must be included and the 
spectra predicted by Monte Carlo must reflect correct 
£'*-decay angular distributions and the effects of ini- 
tial state radiation (ISR). 



4. Evidence for Multi-Body Production 

While the qualitative features of the measured 
charm-meson momentum spectra accorded with ex- 
pectations (Fig. [2), initial attempts to fit the spectra 
did not produce acceptable results. It was quickly con- 
cluded that the two-body processes listed above are 
insufficient to account for all observed charm-meson 
production. Final states like DD^*^mT, in which the 
charm-meson pair is accompanied by one or more 
additional pions, emerged as the likely explanation. 
While not unexpected, these "multi-body" events 
have not previously been observed in the charm- 
threshold region, and there are no predictions of the 
cross sections for Z?" and Z?+ production through 
multi-body final states. 

To assess which multi-body final states {DDtt, 
D*D-K, etc.) are measurably populated in our data 
we examine observables other than the charm-meson 



40 



o 20 



> 



S 40 - 



o 
O 



m 20 



O 



J 

: a) 


n 


ri 

b) 




t — • \ ' \ • 1 — 

- c) 

.7 1.75 1.8 1.85 


J 

1.9 1.95 2 


n : 
J\ - 

2.05 



Missing i\/lass (GeV/c ) 



Figure 3; The mass spectrum of X in (a) e^e D'^-k^X 
at 4170 MeV, (b) e+e" ^ D*^-k'^X at 4170 MeV, and 
(c) e+e" D*°-K^X at 4260 MeV. Peaks at the D* mass 
in (a) and the D mass in (b) are evidence for the decay 
D*Dn. The D peak in (c) confirms D* Dn and the D* 
peak demonstrates that D* D*tv is produced at 4260 MeV. 



momenta, because ISR causes smearing of the peaks 
in the momentum spectra that can obscure the two- 
body kinematics. We applied D^*"^ momentum cuts 
to exclude two-body contributions and examined the 
distributions of missing mass against a and an ac- 
companying charged or neutral pion, using charge cor- 
relations to suppress incorrect combinations. Figure [3] 
shows clear evidence for D*Dtt events at 4170 MeV, as 
well as indications of D*D*tt in the sample of 13 pb~^ 
collected at 4260 MeV (Fig. [3t). These events can- 
not be attributed to two-body production with ISR, 
because radiative photons would destroy any peak- 
ing in the missing-mass spectrum. The absence of a 
peak at the D mass in Fig. [3^ indicates that there is 
no evidence for DDtt production. Analysis of events 
with Ds reveals no evidence for multi-body produc- 
tion, consistent with expectations, since the DfD~Tr^ 
final state violates isospin conservation. 



5. Momentum-Spectrum Fits and 
Cross-Section Results 

Candidate momentum spectra for D'^, D+ and Df 
were selected by requiring the invariant mass to be 
within ±15 MeV of the nominal value. Backgrounds 
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are estimated with a sideband technique. Sideband 
regions are taken on both sides of the expected sig- 
nal, and are significantly larger than the signal region 
to minimize statistical uncertainty in the background 
subtraction. Specific widths are set mode by mode 
based on the expectation of specific background pro- 
cesses. 

Having identified the components of multi-body 
charm production, we determine yields for these chan- 
nels and the two-body modes by fitting the sideband- 
subtracted D'^ and momentum spectra. Sig- 
nal momentum distributions for specific channels are 
based on full GEANT fis'l simulations using EvtGen [l^ 
for the production and decay of charmed mesons. The 
EvtGen simulation incorporates all angular and time- 
dependent correlations by using individual amplitudes 
for each node in the decay chain. ISR is included 
in the simulation, which requires input of energy- 
dependent cross sections for each final state. We used 
simple parameterizations of these cross sections con- 
structed by linearly interpolating between the prelim- 
inary measurements from our analysis. (In doing this 
we make the assumption that the energy dependence 
of the Born-level cross sections is adequately repre- 
sented by the uncorrected cross sections.) For the 
multi-body D*Dtt and D*D*tt final states we used a 
spin-averaged phase-space model within EvtGen. 

Momentum fits for the large sample of data at 
4170 MeV are shown in Fig. H for (a) -> K-n+, 
(b) £)+ K~TT^TT~^, and (c) Df — > 07r+ candidates. 
The lack of entries below ^ 200 MeV confirms 
the absence of multi-body Ds production. Because of 
the relative simplicity of production demonstrated 
by the 4>Tr^ fits and the limited statistics in 

the sample, we determined the final cross sections for 
DfD-, D;+D-, and D;+D;~ by using a sideband 
subtraction technique to count signal events in a re- 
gion of the Afbc and AE plane. The cross sections are 
then determined from a weighted sum of the yields 
for the eight Dg decay modes listed in Sect. [31 with 
weights minimizing the combined statistical and sys- 
tematic uncertainties that were calculated from pre- 
viously measured branching fractions and efficiencies 
determined by Monte Carlo. The cut-and-count anal- 
ysis gives results that are consistent with momentum 
fits. There is good agreement among the separately- 
calculated cross sections for the different Dg decay 
modes. 

Each of the thirteen data subsamples has been an- 
alyzed with the techniques developed and refined on 
data at 4170 MeV. A complete set of fit results is pro- 
vided in Ref. Figure E] shows the D°, D+ and Ds 
fits for data sample at 4260 MeV, which are of par- 
ticular interest because the charm-production cross 
sections might provide insight to the nature of the 
y(4260) state. The fits at 4260 MeV behave similarly 
to those at lower energy, although a larger proportion 
of multi-body decays is apparent. 
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Figure 4: Sideband-subtracted momentum spectra for (a) 
D° ^ K-TV+, (b) D+ K-TV+-K+, and (c) L>+ ^ <;/>7r+ 
(bottom) at 4170 MeV. Data are shown as points with er- 
rors and the total fit result is shown as the solid black line. 
The colored histograms represent fit components, mostly 
single D-production modes. For example, the primary 
in D'^D*'^ , which peaks at 0.7 GeV/c, is shown in bright 
red. The secondary D mesons from the primary D*° de- 
caying via the emission of a tt" form the broad peak at 0.6 
GeV/c shown in light blue. The second broad peak, at 0.6 
MeV/c, consists of D" mesons from the charged pion decay 
of the D*+ in D*'^D~ . All sources of multi-body events 
are combined and result in the broad spectrum between 
and 0.5 GeV/c shown in dark red. 



Cross sections for the two-body and multi-body fi- 
nal states are shown in Fig. [6^-c. The uncertain- 
ties on the data points are statistical and system- 
atic combined in quadrature. Ref. [l7| provides de- 
tailed descriptions of the systematic uncertainties of 
the cross-section determinations. Briefiy, there are 
three sources of systematic uncertainty: determina- 
tion of the efficiency of charm- meson selection, extrac- 
tion of yields, and overall normalization. The total 
systematic uncertainty is not dominated by any one 
of these. 

Track selection and particle identification closely 
follow previous CLEO-c analyses [H, Hi]. The effi- 
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Figure 5: Sideband-subtracted momentum spectra for (a) 
D° K-TV+, (b) D+ K--K+-K+, and (c) Dt 07r+ 
(bottom) at 4260 MeV. Data are shown as points with er- 
rors and the total fit result is shown as the solid black line. 
The colored histograms represent fit components, mostly 
single D-production modes. For example, the primary D° 
in D'^D*^ , which peaks at 0.8 GeV/c, is shown in bright 
red. The secondary D mesons from the primary de- 
caying via the emission of a tt'' form the broad peak at 0.7 
GeV/c shown in light blue. The second broad peak, at 0.7 
MeV/c, consists Z)" mesons from the charged pion decay 
of the D*^ in D*^ D~ . The multi-body events are com- 
bined and result in the broad spectrum between and 0.6 
GeV/c shown in dark red for D* Dn and in black between 
and 0.4 GeV/c for D*D*-k. 



ciency for reconstructing charged tracks has been es- 
timated by a missing-mass technique applied to events 
collected at the ip{2S) and -0(3770) resonances. There 
is good agreement between data and Monte Carlo, 
with an estimated relative uncertainty of ±0.7% per 
track. Pion and kaon identification have been stud- 
ied with and decays in V'(3770) data, with 
estimated systematic uncertainties in the respective 
efficiencies of ±0.3% and ±1.3%. 

The extraction of event yields by fitting the charm- 
meson momentum spectra (non-_Ds modes) incurs sys- 
tematic uncertainty primarily through the signal func- 
tions generated by Monte Carlo, which depend on de- 
tails of ISR and, in the case of D*D*, the helicity am- 
plitudes and resulting _D-meson angular distribu- 
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Figure 6: Exclusive cross sections for two-body and multi- 
body charm-meson final states, and total observed charm 
cross section with combined statistical and systematic un- 
certainties. 



tions. As for the exclusive measurements, these details 
were studied with the large data sample at 4170 MeV, 
for which statistical uncertainties are small, and the 
resulting estimated relative systematic uncertainties 
are applied to all energy points. For the ISR calcu- 
lation, the exclusive cross sections input to EvtGen 
were varied from their nominal shapes. While a qual- 
itative constraint of consistency with our measured 
cross sections was imposed, some extreme variations 
are included in the final systematic uncertainty. Both 
the direct effect on the fitted yield of varying a spe- 
cific mode and the indirect effect of varying other 
modes were computed, although the former dominates 
in quadrature. 

The yields for Dg final states are determined by 
direct counts after cutting on Mbc and AE. Sys- 
tematic uncertainty arises in these measurements if 
the Monte Carlo simulation does not provide an accu- 
rate determination of the associated efficiency. This 
is probed by adjusting the cuts and recomputing the 
cross sections, again using the high-statistics sample 
at 4170 MeV. The systematic uncertainties assigned 
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Table II Total systematic errors in the exclusive cross sec- 
tions. 



Mode 


Relative Error (10 ^) 


Determined by Momentum Fits 




DD 


4.5 


DD* 


3.4 


D*D* 


4.7 


D*Dn 


12.0 




25.0 


Determined by Counting 




DsDs 


5.6 


DsDt 


5.3 




6.8 



based on these studies are ±3%, ±2.5% and ±5% for 
D+D-, D*+D-, and D*+Dl~, respectively. 

In converting the measured yields to cross sec- 
tions we must correct for the branching fractions of 
the charm-meson decay modes. For the non-strange 
charmed mesons, only one mode is used and CLEO- 
c measurements (l^ provide the branching fractions 
and uncertainties: ±3.1% for D'^ K^tt^ and 
±3.9% for D+ K-TT+TT+. For Ds modes we 
use CLEO-c measurements of the branching fractions 
for the eight decay modes included in the weighted 
sum [ijj. The world -average value is used for the 
7?*+ D^TT'^ branching fraction, with a systematic 
uncertainty of ±0.7% Finally, the cross-section 
normalization also depends on the the absolute deter- 
mination of the integrated luminosity for each data 
sample, with a systematic uncertainty of ±1%. 

A mode-by-mode summary of the systematic uncer- 
tainties in the exclusive cross-section measurements is 
provided in Table [TTl 

As a cross-check, for the two largest data samples 
(4170 MeV and 4260 MeV), the multi-body cross sec- 
tions are also determined by fitting the distributions 
of missing mass against detected D^ii, D^tt and D*tt 
combinations. While these measurements are less pre- 
cise, they show good agreement with the results of the 
momentum-spectrum fits. 



6. Inclusive Cross-Section 
Measurements 

If all final states have been included, the sum of the 
exclusive cross sections should equal the total charm 
cross section. We test this supposition with two in- 
clusive measurements that can also be compared with 
past results. 

The first cross-check is a measurement of the total 
charm-meson cross section: 



DDX) 



(1) 



where the contributing cross sections are defined by 
(T£) = N£,leBC, where e and B are the efficiency and 
branching fraction for the decay mode used {D^ — > 
K-t:+, D+ K-t:+tt+, and D+ K-K+TT+), C 
is the integrated luminosity, and Nd is the yield ob- 
tained by fitting the mass spectrum. In the case of 
D*^ and , the invariant-mass distribution is fitted 
to a Gaussian signal and polynomial background. For 
Ds, the event-type requirements are maintained be- 
cause of the relatively large background for the high- 
yield K~K^'!T^ decay mode. For our energy points 
below 4120 MeV, where Ds production occurs only 
through DfDj, the yield is extracted by fitting Aff,c 
to a Gaussian signal and ARGUS background function 
(Tsj . For 4120 MeV and above, event types involving 
D*"*" contribute. For all candidate events that pass the 
selection requirements for any of DfDj , D*^Dj , and 
D;+D;^ (the last only for 4260 MeV), a fit to the Df 
invariant mass is used to determine the yield. 

The second cross-check is a determination of the 
total cross section made by counting multihadronic 
charm events. The contribution of uds continuum 
production is estimated with measurements made at 
Ecm =3671 MeV, below cc threshold, and extrapo- 
lated as 1/s. Procedures for this measurement are 
identical to those used to determine the cross section 
for e+e" ^(3770)^ hadrons in CLEO-c data at 
= 3770 MeV 

Figure shows the inclusive measurements (sta- 
tistical and systematic uncertainties combined in 
quadrature) and the sum of the cross sections for the 
measured exclusive final states, without radiative cor- 
rections. The excellent agreement demonstrates that, 
to current precision, the measured exclusive two- and 
three-body final states saturate charm production in 
this region. Furthermore, charm is demonstrated to 
account for all production of multihadronic events 
above the extrapolated uds cross section. 

For the inclusive-charm cross-section measure- 
ments, the systematic uncertainties associated with 
the per-particle efficiencies for tracking and particle 
identification are identical to those of the exclusive 
measurements. The uncertainties in normalization 
(luminosity and branching fractions) are also identi- 
cal. Systematic uncertainty in the yield extraction is 
dominated by the choice of fitting function. This is 
evaluated mode by mode and propagated into over- 
all systematic uncertainties accounting for all cor- 
relations, with combined systematic uncertainties of 
±4.3%, ±5.1%, and ±8.6% (±10.6%) for D^, D+, and 
Df below (above) 4120 MeV. For the hadron-counting 
inclusive cross sections, the systematic uncertainties 
arc identical to those of Ref . [191 . 

For comparison with other experiments and theory 
it is necessary to obtain Born-level cross sections from 
the observed cross sections by correcting for ISR. We 
do this by calculating correction factors following the 
method of Kuraev and Fadin [So^, which gives the 
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Figure 7: R (with radiative corrections as described in the 
text and correction for non-charm continuum production 
based on Ruds ~ 2.285 ± 0.03, as determined by a ^ fit 
to previous R measurements between 3.2 and 3.72 GeV 
22[) from this analysis and from previous measurements 



observed cross section at any y/s: 
1 

crobs(s) = j dk ■ f{k,s)aB{sca), (2) 



where the Born cross section as is a function of 
the effective center-of-mass energy squared (scff — 
s{l-{E^/E^^^,^))), and f{k, s) is the ISR kernel. The 
radiative-correction factor is also calculated following 
the alternative implementation of Bonneau and Mar- 
tin j2lj . A 4% difference between the two calculations 
is taken as the systematic uncertainty in the radiative 
correction. 

Figure [7] shows that there is excellent agreement 
between our inclusive-charm measurement and R 
measurements in this region made by BES Q and 
Crystal Ball [12]. 



7. Summary and Conclusions 

In summary, we have presented detailed informa- 
tion about charm production above cc threshold. Re- 
alizing the main objective of the CLEO-c scan run, 
we find the center-of-mass energy that maximizes the 
yield of Ds to be 4170 MeV, where the cross section 
of ~ 0.9 nb is dominantly D*^D~ . This informa- 
tion has guided the planning of subsequent CLEO- 
c runnin g, w ith initial results already presented on 
leptonic ]26| and hadronic Dg decays. The to- 
tal charm cross section between 3.97 GeV and 4.26 
GeV has been measured both inclusively and for spe- 
cific two-body and multi-body final states. Internal 
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Figure 8: Comparisons between measured cross sections 
and the updated predictions of the potential model of 
Eichten et al. [ii2^] (solid lines). 

consistency is excellent and radiatively-corrected in- 
clusive cross sections are consistent with previous ex- 
perimental results. Figure [6] shows that the observed 
exclusive cross sections for DD, D*D, D*D* , DfD~, 
DI+D-, D*+D*-, D*Dtt, and D*D*t: exhibit struc- 
ture that reflects the intricate behavior expected in 
the charm-threshold region. Figure [5] provides a com- 
parison between our measured cross sections and the 
updated calculation of Eichten et al. [2^ . There is 
reasonable qualitative agreement for most of the two- 
charm-meson final states. The most notable exception 
is the cross section for D*D* in the region between 
4050 and 4200 MeV, where the measurement exceeds 
the prediction by as much as 2 nb. This corresponds 
to nearly a factor-of-two disagreement in the ratio of 
D*D* to D*D production, accounting for about two 
thirds of the difference in the total charm cross sec- 
tion. This is a much larger effect than the absence of 
a multi-body component from the theoretical predic- 
tion. 

It has been suggested by Dubynskiy and Voloshin 
[25} that the existence of a peak in the D*D and DsDg 
channels at the D*D* threshold, along with the obser- 
vation that there is a minimum in DD, can be inter- 
preted as a possible new narrow resonance, but avail- 



8 



Proceedings of the CHARM 2007 Workshop, Ithaca, NY, August 5-8, 2007 



able data are insufficient for a definitive assessment. 

The D*D* cross section exhibits a plateau just 
above its threshold. This contrasts with D*D, which 
we observe to peak at threshold, in agreement with 
recently presented preliminary results from Belle (23j . 
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